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Abstract. We report on new data, including four vacuum ultravio- 
let spectropolarimetric observations by the Wisconsin Ultraviolet Photo- 
Polarimeter Experiment ("WUPPE") on the Astro-1 and Astro-2 shut- 
tle missions, and 15 new visible-wavelength observations obtained by the 
HPOL CCD spectropolarimeter at the Pine Bluff Observatory of the Uni- 
versity of Wisconsin. This includes three HPOL observations made within 
12 hours of each of the three Astro-2 WUPPE observations, giving es- 
sentially simultaneous observations extending from 1500 to 10500A An 
analysis of these data yields estimates on the properties of wind "clumps" 
when they are detectable by the polarization of their scattered light. 
We find that the clumps must be detected near the base of the wind, 
r/i?* = 1.3 — 2.5. At this time the clump density is at least 10 13 cm -3 , 20 
times the mean wind, and the temperature is roughly 10,000 °K, about 
20% cooler than the mean wind. The clumps leading to the largest ob- 
served polarization must have electron optical depth of 0.1 - 1 and an 
angular extent of 0.1 - 1 ster, and account for at most 2% of the wind 
mass loss. The H and Hel emission lines from the wind are unpolar- 
ized, but their P Cygni absorption is enhanced by a factor of four in the 
scattered light. We speculate on the possible relationship of the clumps 
detected polarimetrically to those seen as Discrete Absorption Compo- 
nents (DACs) and in interferometry. 



1. Introduction 



P Cygni was found by many early polarimetric investigators to possess intrinsic 
linear polarization (eg, Coyne & Gehrels 1967; Coyne, Gehrels, & Serkowski 
1974; Serkowski, Mathewson, & Ford 1975). Hayes (1985) and Lupie & Nord- 
sieck (1987) established that the polarization was variable on timescales of days 
- weeks, with an amplitude of approximately 0.4%, and with no apparent peri- 
odicity nor favored position angle, although the interpretation was hampered by 
an appreciable, unknown interstellar polarization. The long time series of Hayes 
(1985) exhibits a polarimetric episode length of ~ 1 week, and a repetition rate 
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of once per 20-30 days. Taylor, Nordsieck et al (1991, hereafter "TN"), in the 
first detailed spectropolarimetric investigation, made use of the apparent non- 
variability of the polarization of the Ha emission line to estimate the interstellar 
polarization there, extrapolating to other wavelengths using a galactic mean in- 
terstellar wavelength dependence. After removing this interstellar polarization, 
TN then confirmed the remarkable daily variability and the random nature of 
the degree and position angle of the intrinsic continuum polarization. The Ha, 
HP, and Hel A5876 lines were found to be essentially unpolarized, while the 
continuum polarization wavelength dependence is not flat, but exhibits a steady 
decrease from the near ultraviolet into the red. TN pointed out that this, cou- 
pled with the lack of a polarization Balmer Jump, was difficult to explain using 
the model which is conventionally used for Be stars, where hydrogen bound- free 
opacity competes with wavelength-independent electron scattering to impart a 
polarization wavelength dependence. 

2. New Spectropolarimetry 

The new polarimetric data presented here includes four vacuum ultraviolet spec- 
tropolarimetric observations by WUPPE on the shuttle missions Astro-1 (one 
observation: 12 Dec, 1990) and Astro-2 (three observations: 3, 8, and 12 March, 
1995), and 15 new visible- wavelength observations obtained by the "HPOL" 
CCD spectropolarimeter at the Pine Bluff Observatory of the University of 
Wisconsin. The WUPPE instrument is described in Nordsieck et al (1994) 
and HPOL is described in Nordsieck & Harris (1995). The visible wavelength 
spectropolarimetry includes three HPOL observations made within 12 hours of 
each of the three Astro-2 WUPPE observations, giving essentially simultaneous 
observations extending from 1500 to 105004 . The Astro-1 WUPPE observation 
was discussed by Taylor, Code, et al (1991), who noted a total continuum polar- 
ization generally consistent with interstellar, plus possible line-blanketing of the 
intrinsic polarization, an effect similar to that observed in WUPPE observations 
of Be stars (Bjorkman, et al 1991). 

Figure 1 shows the net polarization vs wavelength from the WUPPE obser- 
vations and the three contemporaneous HPOL observations. The observations 
exhibit an overall interstellar polarization continuum shape plus strong intrinsic 
polarization variability over most of the spectrum. The interstellar polarization 
is estimated in a variety of ways. First, we see that the four WUPPE observa- 
tions all cross in the 1700 — 19004 region. We interpret this to be strong Felll 
line blanketing which is suppressing the intrinsic polarization there. This region 
is known to be the most strongly blanketed region in P Cyg (Pauldrach and Puis 
1990), and strong polarization suppression there is seen similarly in the Be star £ 
Tau (Bjorkman, et al 1991). Conveniently, this then establishes the interstellar 
polarization value at 18004, 0.47 ± 0.03% at PA 32.1. Second, following TN, 
the Ha line is assumed to be not intrinsically polarized. This is supported by 
figure 2a, which shows a polarization vector diagram for all of the HPOL ob- 
servations. The continuum polarization near Ha varies by roughly ±0.3% with 
no apparent favored position angle, while the Ha points are consistent with a 
constant polarization of 1.11 ±0.03% at PA 35.5. This is then taken to be the in- 
terstellar polarization at Ha (cross in figure 2a). Using the standard Serkowski 
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Figure 1. Observed polarization of P Cyg. Top: flux. Middle: degree 
of linear polarization, for the Astro-1 and the three Astro-2 observa- 
tions, with corresponding HPOL observations. Bottom: polarization 
position angle. Gray dashed line: interstellar polarization estimate 
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Figure 2. 2a (left): Polarization vector diagram for 15 HPOL obser- 
vations. 2b(right): Interstellar polarization map 
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Figure 3. P Cyg mean intrinsic polarization wavelength depen- 
dence. Fit with competing free-free absorption for clump optical depth 
Tff(lfJLm) = 1, 1.5, and 2. 

interstellar polarization law as modified by Wilking, Lebowski, & Rieke (1987; 
"WLR"), plus a modest position angle rotation, these two points determine the 
actual P Cyg WLR parameters: p max = 1.17 ± 0.03%, X max = 5100 ± 100A, 
P. A. = 35.1°, and a position angle rotation of —0.8 ± 0.5°//mi _1 . The adopted 
interstellar polarization curve is shown as the gray dashed line in figure 1. We 
have attempted to verify the above estimate through observations of stars near 
P Cyg (figure 2b). Observations with HPOL with the 0.9m telescope at PBO of 
HD19300 and HD19307 (17 arcmin distant) exhibit rather different interstellar 
polarizations, consistent with the large known variations in Cygnus. New obser- 
vations at the WIYN 3.5m telescope of two B stars less than 5 arcmin distant 
within the P Cyg cluster, LSII 37-49 and -54, yield more comparable results, 
1.5% at 35° and 0.9% at 47°. This provides support for our adopted value. 



3. Intrinsic Polarization 

In Figure 3 we show the mean intrinsic polarization wavelength dependence 
for the three Astro-2 + HPOL observations. This was obtained by removing 
the fitted interstellar polarization, rotating the intrinsic polarization into the Q 
stokes parameter using the blue position angle, and averaging the result. This 
is the same procedure used by TN to estimate the mean wavelength dependence 
of the intrinsic polarization: it assumes that the intrinsic polarization at one 
time has a position angle independent of wavelength. The Stokes polarization Q 
is positive when the polarization is parallel to the blue continuum polarization, 
and negative when it is orthogonal. Note the strong decrease into the vacuum 
ultraviolet (presumably due to line blanketing) , the strong decrease into the red 
as noted by TN but here extended into the near infrared, and the lack of both 
a B aimer and Paschen jump . 
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3.1. Clump Optical Depth 

It is assumed that for a star as hot as P Cyg, any continuum intrinsic polar- 
ization is due to electron scattering in the envelope, presumably in asymmetric 
inhomogeneities ("clumps"). We may estimate the electron scattering optical 
depth r e of the clumps from a simple single-scattering analysis of the degree of 
polarization. Using formulae in Cassinelli, Nordsieck, and Murison (1987), the 
net polarization from a sector- shaped density enhancement subtending An < 1 
ster and with electron scattering optical depth r e < 1 is 

p » (47r) _1 /( 1 + (47T)" 1 ) sin 2 i D(r) r e An 
= 0.073 sin 2 i (1 - (R pho t/r) 2 ) 1/2 Te An. 

Here i is the inclination of the clump direction to the line of sight, and D(r) is 
the "finite disk correction factor" which corrects for cancellation of polarization 
near the stellar photosphere caused by illumination from a non-point source 
(here we take the approximation of no photospheric limb darkening). Given 
the rapid timescale and the large amount of polarization observed, we expect 
that the clumps are polarimetrically observable only near the base of the wind. 
Later we find that the appropriate value for the radial position of the clumps is 
r/R* = 1.4 - 2.5, which gives D ~ 0.4 - 0.9 for R pho t/R* = 1-3 (Najarro 2000); 
we will take D = 0.65. This formula holds approximately until r e ~ 1, where 
the clump becomes optically thick to the total opacity, and An ~ 1 ster, where 
geometric cancellation over the clump begins to reduce the polarization. Thus 
the maximum possible polarization (for i = 90) is 

Pmax ~ 0.073 x 0.65 r e Afi ~ 0.05 r e An. 

The observed polarization amplitude below 5000A is roughly 0.004, so that 
r e An ~ 0.1. A range of clump properties consistent with the observations is 
then (r e , An) = (0.1, 1) to (1, 0.1). 

3.2. Free-Free Absorption - Clump Density 

The presence of continuum absorption optical depth r a in addition to the electron 
scattering optical depth r e reduces the net polarization by the absorption in the 
clump 

p(X) = F p (X)/F(X) ex r e e~ Ta 

where F is the total flux and F p is the polarized scattered flux. Features in 
p(X) may then be used to estimate r (A), and thus the physical properties of the 
clumps. 

An important feature of the observed continuum intrinsic polarization wave- 
length dependence is the persistent decrease into the infrared. This was also 
noted by TN, and holds for individual nights as well as the mean. A plausible 
explanation for this, given the high density and ionization at the base of the 
wind, is competing free-free absorptive opacity. Since kjj oc A 2 , the polarization 

will decrease into the red as p(X) oc e~ aA2 . Figure 3 shows a fit to the polariza- 
tion for Tff(l(j,m) = 1, 1.5, and 2. We adopt 1.5 as the best value. Given our the 



6 



Nordsieck et al 



range of allowed r e = 0.1 — 1, we have rgr(l/xm)/r e = 1.5 — 15. This provides an 
estimate of the clump electron density: 

T ff (X)/T e = K ff (X)/a = a ff X 2 gN e N p T- 3 / 2 /N e a T (cgs units), 
or N e (clump) = 3 x 10 12 T 4 (clump) 3/2 T ff {l fim)/r e cm~ 5 ~ 10 13 cm" 3 . 

(For these rough estimates we assume a highly ionized pure hydrogen gas, T4 is 
the clump temperature in units of 10 4o K, and g ~ 1 is the Gaunt factor). This 
is approximately 20 times the density in the wind at these radii in the model of 
Najarro (2000). 

3.3. Bound-Free Jumps: Clump Temperature and Ionization 

The lack of Balmer or Paschen Jumps in the intrinsic polarization of figure 
3 puts an upper limit T6/(Ab = 3646) < 0.1 for neutral hydrogen bound-free 
absorption in the clumps, which then means T6y(As)/rgr(l/im) < 1/15. This 
ultimately puts a limit on the clump temperature: 

n f (X B ) _ a hf (X B )N H (n = 2) _ a bf (X B ) e ~W-2/kT 
T ff (lnm) a ff (lnm)gN e N p T-V 2 a ff (lfim) H 2 N e T~ 3 / 2 ' 

where Xh is the neutral fraction and 62 is the departure coefficient for the 
second state of hydrogen in the clump. For simplicity we will assume that 
62 is unchanged in the clump; it is less than two at the base of the wind in 
the models of Najarro (2000), in any case. For Xh, we assume pure optically 
thin photoionization, so that it may be scaled from the base wind density and 
temperature: 

N(wind) \ T(wmd) ) 

Finally, then, 

n/ (Ab) =75x 1Q i9 b 2 X H (wind) T 4 {wind) - 8 ^ 
Tff(lfj,m) N(wind) 

T 4 (clump)°- 7 e- n - 8 / T ^ dum ^ < 1/15 

Given a wind model, this may be solved for T (clump), giving the maximum 
temperature of the clump consistent with no detectable bound-free absorption. 
Notice that this is independent of clump electron scattering optical depth and 
density. Figure 4 (heavy lines) shows the clump density and ionization corre- 
sponding to the maximum temperature at each radius in the wind, using the 
wind model of Najarro (2000), for r e = 0.1 and 1. The clumps will not be de- 
tectable inside r = 1.312*, the position of the photosphere. Above r = 2.512*, 
the maximum temperature decreases below 8, 000° K and recombination must 
occur; this puts a rough upper limit on the radius. 

The clumps must contribute their maximum polarization from r = 1.3 to 
2.512*, and at this time they must be at least 20 times more dense and 20% 
cooler than the wind. 
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Figure 4. Possible polarimetric clump properties as a function of 
radial position. Left: density (solid lines, lefthand scale); temperature 
(dashed lines, righthand scale). Right: neutral fraction (solid, lefthand 
scale); wind velocity (dashed line, righthand scale). For both, wind: 
light lines; clumps: heavy gray lines (r e = 0.1), heavy black lines (r e = 
1). The clump temperature is the maximum allowed for all clump 
optical depths. 



3.4. Clump Geometry and Time Evolution 

It is interesting to compare the geometric properties of the range of possible 
clumps. "Thin" clumps have r e = 0.1 and AO = 1. "Thick" ones have 
r e = 1 and AO = 0.1. The physical thickness of these clumps is Ar = 
T e/ \&TN e (clump)) = 0.0005 and 0.05i2*, respectively, at r/R* = 1.5. Figure 5 
illustrates this " pancake- like" geometry. The mass of the thick and thin clumps 
is the same, a result of the requirement that they produce the same polarization. 
The instantaneous stellar mass loss rate due to a clump compared to that in the 
base wind is (JV (clump) /N (wind)) AQ / 4 vr = 0.16 and 16 at r/R* = 1.5, which 
lasts for a time Ar/v = 1 and 0.01 day, respectively. Given that the repetition 
rate for polarimetric clumps is 20 - 30 days, these clumps account for at most 
1% of the mass loss of the star. The largest mass clumps are required at the 
largest distance from the star: at r/R* = 2.5, they still account for at most 2% 
of the wind mass loss. 

Finally, we may speculate on the time evolution of the clumps and their 
polarization. This will depend on the evolution of the mass and density of the 
clumps as they are carried out in the wind. The analysis above has placed 
limits on the properties of the clumps at the time that they contribute the 
largest polarization. The polarization will grow as the clump appears out of 
the photosphere, possibly growing in mass as it forms. It may then dissipate 
and its density, optical depth, and polarization will drop as it is carried out in 
the wind. This is a complex problem requiring a proper radiative transfer and 
clump physics model. We may explore the purely geometric effects, however: 



8 



Nordsieck et al 




phqto'sphere 



Figure 5. Clump cartoon, to scale. 

if we ignore optical depth effects and changes in clump mass, p(r) oc D(r)/r 2 . 
This has a maximum at r/R p h t ~ 1-1 (r/R* ~ 1-4), and drops by a factor of 10 
by r/i?# = 3. This is consistent with the polarizing clumps being highly ionized: 
by the time they reach the recombined part of the wind, the polarization is no 
longer detectable. The time it takes for a clump to traverse the r/R* = 1.3 — 2.5 
region will be 1.2i2*/(55km/s) = 14 days, which compares favorably with the 
observed duration of roughly one week. 

3.5. Intrinsic Line Polarization 

Figure 6 shows an expanded version of the visible-wavelength mean intrinsic 
polarization spectrum from all 15 HPOL observations. The improved sensitivity 
of the new HPOL data allows one to extend the TN analysis to 4 Hydrogen 
lines and 5 Hel lines. The intrinsic polarization of the emission line flux itself 
is a diagnostic of the wind. Any line radiation produced deep in the wind may 
become polarized if it is scattered off the clumps that produce the continuum 
polarization. Alternatively, the P-Cygni absorption beneath the emission line 
(unresolved in these observations) may be different in the clump, producing an 
apparent line polarization. The latter appears to be the largest effect in P Cyg. 
The line Stokes polarization qi is conventionally computed as follows: 

qi = J (J(A)g(A) - I c q c ) dX I J (/(A) - I c ) dX 

where the integral is performed over the observed line profile, and I c and q c are 
the flux and Stokes polarization in the neighboring continuum. The initially 
surprising result is that the weaker lines all appear to be negatively polarized 
(polarized orthogonal to the continuum). We believe that this is an artifact of 
the fact that the P Cygni absorption is enhanced in the clumps, and is as a 
result interpreted as emission polarized orthogonal to the continuum. For such 
a situation one can show that the apparent net Stokes polarization qi of the line 
is 

_ qe-fqcWg/We 
qi l-W a /W e ' 

where q e and q c are the Stokes polarization of the emission line and continuum, 
respectively, W a and W e are the equivalent widths of the absorption and emission 
part of the P-Cygni line, respectively, and / = W a (pol)/W a is the factor by 
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Figure 6. P Cyg visible wavelength mean polarization wavelength 
dependence, showing polarimetric line effects and line identifications. 
Top: flux. Second panel: degree of polarization. Third panel: position 
angle. Bottom: "polarized flux", showing enhanced absorption. 



which the absorption line is enhanced in the polarized light. W a and W e were 
evaluated using the mean P Cyg line profiles in Stahl, et al (1993). 

In Figure 7 we plot qi (1— W a /W e ) versus q c W a /W e for the 9 lines measured. 
In this plot, the intercept of a straight line fit to the lines of one species is q e and 
the slope is /. The line labeled "unpol" in this diagram corresponds to q e = 
and / = 1. We see that both the H and Hel lines are consistent with unpolarized 
emission (q e = ± 0.02%), implying that more than 90% of the emission arises 
outside of the polarizing region. On the other hand, both H and Hel are strongly 
enhanced in absorption in the polarized flux (f ~ A). The enhancement must 
be due to a complex combination of ionization and excitation effects, and will 
require a detailed analysis. It is perhaps surprising that the values are so similar 
for H and Hel. Another way of seeing these effects is in the bottom panel of figure 
6, the "polarized flux". This is essentially the spectrum of the scattered light: 
the emission lines have almost vanished, indicating the emission is unscattered, 
and strong absorption is seen instead. 



4. Relation of Polarimetry to other Observations 

A frequently mentioned datum for P Cyg is the apparent rotation from vsvai 
measurements. Markova (2000) favors a value of 40 km/s, and it is agreed that 
values less than 25 km/sec are not consistent with absorption line depths. 40 
km/s is 15% of breakup velocity for a 30Mq, 75-Rq star (270 km/s). This 
is large enough that one might expect a substantially flattened wind if viewed 
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Figure 7. Line polarization diagnostic diagram. Solid diamonds: H; 
Crosses Hel (ortho); squares: Hel (para). Dashed line: locus for unpo- 
larized emission and unenhanced absorption. Solid and dashed lines: 
unpolarized emission plus enhanced absorption for H and Hel. 



edge-on. However, the intrinsic polarization does not show a favored position 
angle: In figure 2a the continuum polarization points occupy a nearly circular 
region. By contrast, Be stars, which are rotating at ~ 50% break-up, show 
polarization vector diagrams confined to a narrow line. One way out of this 
apparent contradiction is if P Cyg is rotating rapidly, and is viewed not quite 
pole-on. For such an object, assuming that all clumps are ejected at random in 
the equatorial plane, the expected axial ratio in the polarization vector diagram 
is ~ cos 2 i. If we insist that cos 2 i > 0.9, this gives i < 18°. This in turn gives 
v > wsini/0.3 = 120 km/s, which is still 45% breakup. This would likely be a 
bipolar object, but viewed nearly pole-on. Several things recommend this sce- 
nario: First, P Cyg has an unusually high wind speed among S Dor variables 
(Luminous Blue Variables). This would be explained if we are just observing 
a fast polar wind in a bipolar object. Second, most S Dor's are bipolar. And 
third, most clumps would be ejected near the plane of the sky, which is favorable 
for large polarimetric variations. On the other hand, spectroscopy of the P Cyg 
shell (Meaburn 2000) does not show the characteristic multiple- redshift compo- 
nents seen in bipolar Planetary Nebulae, for instance. The pole-on hypothesis 
is testable by looking for a systematic daily change in the polarimetric position 
angle. During the 14 days of detectability, a clump would rotate about 60° about 
the star, and each successive clump would exhibit the same sense of rotation. 

How might the polarimetric clumps be related to the Discrete Absorption 
Components (DAC's) seen in P Cyg? If there were no rotation, or the object is 
viewed pole-on, they would probably be uncorrelated, since polarimetric clumps 
must be ejected near the plane of the sky and would never cover the stellar disk 
to form DAC's. On the other hand, if the star is seen equator-on and is rotating 
at v ro t = 40 km/s, a clump will have rotated about 20°by the time it reaches 2i?*, 
and 60° by 10-R* (if the clump maintains its integrity this long). This suggests 
that there would be a weeks - month time delay between the appearance of a 
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polarimetric clump and its observability as a DAC. This correlation would vary 
with the initial ejection longitude of the clump. Given the relatively frequent 
generation of polarimetric clumps (20-30 days), such a variable delay might be 
difficult to establish. 

On the other hand, the polarimetric clumps should certainly be related to 
Ha clumps observed interferometrically (Chesneau, et al 2000). If the clump 
does not disintegrate first, there should be a roughly 2-6 week delay between 
the appearance of a polarimetric clump and its detectability as an Ha clump 
at 10-R*. The position angles of the polarimetric and interferometric clumps 
should be correlated, and systematic differences in position angle would be a 
sensitive measure of the component of stellar rotation about the line of sight 
(which combined with v sin i would yield the inclination) . 

5. Summary; Future 

In summary, we have analyzed the variability of the wavelength dependence of 
the polarization of P Cyg from 1500 — 10500A . After removal of the interstellar 
polarization, we have put limits on the properties of the wind inhomogeneities 
("clumps") during the brief period that they are detectable in polarization. 

• The presence of free-free absorption in the polarization wavelength de- 
pendence implies a clump electron density of at least 10 13 cm -3 . This 
corresponds to a density contrast of at least 20 at the base of the wind. 

• The absence of Balmer and Paschen bound-free absorption in the polar- 
ization limits the polarimetric zone to 1.3 — 2.5i?*, where the hydrogen 
ionization is high, and requires the clumps to be substantially cooler than 
the wind. The clumps take 14 days to be swept through this region. 

• In order to produce the maximum observed polarization, the largest clumps 
must have an electron optical depth of 0.1 - 1 and an angular extent of 0.1 
- 1 ster, accounting for at most 2% of mean mass loss of the wind. 

• The emission lines are essentially unpolarized. However, the P Cygni ab- 
sorption is enhanced by a factor of four in the polarimetric clumps. 

We plan on the following further observational investigations: 

• Higher spectral resolution R ~ 3000 spectropolarimetry to resolve polari- 
metric line profiles. This should verify the line analysis model of section 
3.5. 

• Monitor the polarization every night for at least a month. During such a 
campaign we would look for a change in polarization wavelength depen- 
dence during an episode, including the expected decease in the free-free 
absorption and the possible appearance of bound-free jumps near the be- 
ginning or end of the episode. We would also look for a systematic position 
angle change with time, indicating rotation. We will attempt to arrange 
simultaneous high-resolution spectroscopy (for DACs) and interferometry, 
to pursue the predictions of section 4. 
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For modeling, we plan to develop our Monte Carlo polarization radiation 
transfer code to give a quantitative continuum polarization model based on the 
mean wind from Najarro (2000), including electron scattering, H bound free, H 
free-free, and line blanketing, with a simple sector geometry. This would predict 
the absolute degree of polarization, the continuum wavelength dependence, and 
their time dependence as the clump is carried out in the wind. Variables would 
include the clump heating/ cooling and confinement mechanisms. 

A more ambitious polarization model would treat the lines. To accomplish 
this, it will be necessary to more correctly model the ionization and populations 
in the clump. Such a model would predict equivalent width in polarized flux for 
H and Hel to compare with the observations. Finally, our ultimate goal is to cal- 
culate polarimetric line profiles, allowing for electron scattering broadening and 
Doppler shifts in the wind, to compare with future high resolution observations. 
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